
import numpy as np
import matplotlib.pyplot as plt

p53seg = "cggagcagctcactattcacccgatgagaggggaggagagagagagaaaatgtcctttaggccggttcctcttacttggcagagggaggctgctattctccgcctgcatttctttttctggattacttagttatggcctttgcaaaggcaggggtatttgttttgatgcaaacctcaatccctccccttctttgaatggtgtgccccaccccccgggtcgcctgcaacctaggcggacgctaccatggcgtagacagggagggaaagaagtgtgcagaaggcaagcccggaggcactttcaagaatgagcatatctcatcttcccggagaaaaaaaaaaaagaatggtacgtctgagaatgaaattttgaaagagtgcaatgatgggtcgtttgataatttgtcgggaaaaacaatctacctgttatctagctttgggctaggccattccagttccagacgcaggctgaacgtcgtgaagcggaaggggcgggcccgcaggcgtccgtgtggtcctccgtgcagccctcggcccgagccggttcttcctggtaggaggcggaactcgaattcatttctcccgctgccccatctcttagctcgcggttgtttcattccgcagtttcttcccatgcacctgccgcgtaccggccactttgtgccgtacttacgtcatctttttcctaaatcgaggtggcatttacacacagcgccagtgcacacagcaagtgcacaggaagatgagttttggcccctaaccgctccgtgatgcctaccaagtcacagacccttttcatcgtcccagaaacgtttcatcacgtctcttcccagtcgattcccgaccccacctttattttgatctccataaccattttgcctgttggagaacttcatatagaatggaatcaggatgggcgctgtggctcacgcctgcactttggctcacgcctgcactttgggaggccgaggcgggcggattacttgaggataggagttccagaccagcgtggccaacgtggtg"

standard = { 'ttt': 'F', 'tct': 'S', 'tat': 'Y', 'tgt': 'C',
		'ttc': 'F', 'tcc': 'S', 'tac': 'Y', 'tgc': 'C',
		'tta': 'L', 'tca': 'S', 'taa': '*', 'tga': '*',
		'ttg': 'L', 'tcg': 'S', 'tag': '*', 'tgg': 'W',

		'ctt': 'L', 'cct': 'P', 'cat': 'H', 'cgt': 'R',
		'ctc': 'L', 'ccc': 'P', 'cac': 'H', 'cgc': 'R',
		'cta': 'L', 'cca': 'P', 'caa': 'Q', 'cga': 'R',
		'ctg': 'L', 'ccg': 'P', 'cag': 'Q', 'cgg': 'R',

 		'att': 'I', 'act': 'T', 'aat': 'N', 'agt': 'S',
 		'atc': 'I', 'acc': 'T', 'aac': 'N', 'agc': 'S',
		'ata': 'I', 'aca': 'T', 'aaa': 'K', 'aga': 'R',
  		'atg': 'M', 'acg': 'T', 'aag': 'K', 'agg': 'R',

		'gtt': 'V', 'gct': 'A', 'gat': 'D', 'ggt': 'G',
		'gtc': 'V', 'gcc': 'A', 'gac': 'D', 'ggc': 'G',
		'gta': 'V', 'gca': 'A', 'gaa': 'E', 'gga': 'G',
		'gtg': 'V', 'gcg': 'A', 'gag': 'E', 'ggg': 'G'
		}

length = len(p53seg)
i =0
gccount=0
while i<length:
	if (p53seg[i] == 'g') | (p53seg[i]=='c'):
		 gccount+=1
	
	i+=1
i-=1
ratio=1.0*gccount/len(p53seg)
print(ratio)
reverse = p53seg
reverse = reverse.replace('a', 'b')
reverse = reverse.replace('g','d')
reverse = reverse.replace('t', 'a')
reverse = reverse.replace('c', 'g')
reverse = reverse.replace('b','t')
reverse = reverse.replace('d', 'c')
reverse = reverse[::-1]
print(reverse)
print(p53seg[1])
print(standard['gta'] + standard['gcc'])

print("Frame One")
j=0
codon = 'aaa'
amino = " "
while j < len(p53seg)-2:
	codon = p53seg[j] + p53seg[j+1] + p53seg[j+2]
	amino = amino+ standard[codon] 
	j+=3

print(amino)

